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SIMPLE OPTIMAL ESTIMATOR FOR 
PbA STATE OF CHARGE 

FIELD OF THE INVENTION 
[0001] The present invention relates to power sources, and more 
particularly to a state of charge estimator for a lead-acid (PbA) power 
source. 

5 

BACKGROUND OF THE INVENTION 
[0002] The state of charge (SOC) of a battery indicates the amount 
of usable energy stored within the battery at a given time. It is 
analogous to the amount of fuel in a fuel tank. In order to improve 

10 battery life and overall battery performance, the battery must be 

maintained at an optimal SOC. The state of health (SOH) of a battery 
indicates the condition of the battery. The SOH is based on observing 
particular battery parameters to determine the extent of unobservable, 
internal damage or degradation. 

15 [0003] A battery can be considered a system that has current and 
heat flow as inputs and responses of voltage and temperature. The 
voltage response of the battery to current stimulation is described using 
a combination of voltages and a current source. These voltages 
include an equilibrium potential, voltage resulting from hysteresis 

20 effects, voltage resulting from ohmic resistance drops, and voltage 

resulting from dynamics in the battery, or double-layer voltage. Each of 
these voltages can be described by either an algebraic function, a 
differential equation or a convolution integral. The current source is 
indicative of a self-discharge current of the battery. 

25 [0004] During use, battery resistance, equilibrium potential, voltage 
resulting from hysteresis effects, voltage resulting from ohmic 
resistance drops, double-layer voltage, and self-discharge current are 
not directly measurable. These values influence the SOC and the 
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parameters that indicate the SOH of the battery. Because they are not 
directly measurable, it is often difficult to accurately determine the SOC 
and SOH of the battery. Therefore, it is difficult to maintain the battery 
at an optimal SOC or to determine when the SOH has degraded to a 
5 point where battery performance is significantly effected. 

SUMMARY OF THE INVENTION 
[0005] Accordingly, the present invention provides a state of charge 
estimator for a battery. The state of charge estimator includes a meter 

10 that generates a terminal voltage signal of the battery and a terminal 
current signal of the battery. A controller employs a linearized model of 
the battery and a time-varying state estimator. The controller 
processes a synthesized input based on the terminal current and the 
terminal voltage to estimate the battery state of charge. 

15 [0006] In one feature, the battery system includes a generator that 
is controlled by the controller to selectively charge the battery based on 
the state of charge. 

[0007] In another feature, the controller predicts a current state of 
charge based on a prior state of charge and predicts a current state of 
20 charge error based on a prior state of charge error. The controller 
determines a current gain based on the current state of charge error 
and updates the current state of charge based on the prior state of 
charge. 

[0008] In another feature, the synthesized input compensates for 
25 gassing and self-discharging of the battery. 

[0009] In still another feature, the synthesized input compensates 

for dynamic impedances and non-linear characteristics of the battery. 

[0010] In yet another feature, the battery is a lead-acid battery. 

[001 1 ] Further areas of applicability of the present invention will 
30 become apparent from the detailed description provided hereinafter. It 

should be understood that the detailed description and specific 

examples, while indicating the preferred embodiment of the invention, 
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are intended for purposes of illustration only and are not intended to 
limit the scope of the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 
5 [0012] The present invention will become more fully understood 
from the detailed description and the accompanying drawings, wherein: 
[0013] Figure 1 is a functional block diagram of a system including a 
battery, a load powered by the battery and a generator to recharge the 
battery; 

10 [0014] Figure 2 is a schematic diagram of the battery; 

[0015] Figure 3 is a graph of state of charge (SOC) based on open 
circuit voltage (V Q ); 

[0016] Figure 4 is a flowchart illustrating the state of charge 
estimator control according to the present invention; and 
15 [0017] Figure 5 is a flowchart illustrating estimation control for the 
state of charge estimator control of Figure 4. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[0018] The following description of the preferred embodiment is 
20 merely exemplary in nature and is in no way intended to limit the 
invention, its application, or uses. For purposes of clarity, the same 
reference numbers will be used in the drawings to identify similar 
elements. 

[0019] Referring now to Figure 1 , a system 10 includes a battery 12 
25 that is discharged as it powers a load 14. The battery 12 is selectively 
recharged by a generator 16. Operation of the generator 16 is 
controlled by a controller 18. A meter 20 communicates with the 
battery 12 and generates a terminal voltage signal (V t ) and a terminal 
current signal (l t ). The controller 18 process the terminal voltage and 
30 terminal current signals according to the present invention to estimate 
a state of charge (SOC) of the battery 12. The controller 18 signals the 
generator 16 to operate based on the SOC estimate. 
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[0020] Although the system 10 is a generic battery system, it is 
indicative of a battery system implemented in a vehicle. In a vehicle, 
the battery 12 powers loads such as headlamps, heaters, power seats 
and windows, and the like. The vehicle's engine drives the generator 
5 16 or alternator to produce recharging current. 

[0021] Referring now to Figure 2, the battery 12 is schematically 
modeled in terms of an equilibrium potential (V Q ), voltage resulting from 
hysteresis effects (V hys t), voltage resulting from ohmic resistance drops 
(Vohmic), voltage resulting from dynamics in the battery (V d i), or double- 

10 layer voltage and, and a self-discharge current (I S d). V q can be 

determined based on V dl and Vhyst- Vdi and Vhyst sre determined using 
an open loop model. These voltages are subtracted from a 
measurable terminal voltage (V t ) to determine V G . Knowing V Q is 
sufficient to estimate the SOC of the battery 12 given a function that 

1 5 maps V 0 to SOC. An empirical fit of the SOC to V Q for an exemplary 
lead-acid (PbA) battery is shown in Figure 3. Although the battery 
model of Figure 2 does not include all of the existing battery behaviors, 
testing and experience has shown that this battery model is sufficient 
for SOC estimation, power control and diagnostics for automotive 

20 applications. 

[0022] The present invention provides a Kalman Filter (KF)-based 
state estimator to determine the SOC and resistance of the battery 12. 
A standard Kalman Filter (KF) combines the concept of recursive least 
squares estimation with a state space model and noise model to 

25 optimally estimate the state (x) of a system. The KF estimates the 

state of the system based on a linear model of the system. The linear 
model of the system includes a linear dynamic model, a description of 
the covariance of the errors in the system dynamics and a description 
of the covariance of the errors in system measurements. The linear 

30 model is expressed in state space form as: 

x k+ i = A k x k + B k u k + w k (Process Model); and 
y k = F k x k + G k u k + r k (Measurement Model) 



4 



where: 

x k is the state vector (i.e. what is to be estimated); 
y k is the measurement vector (known); 
u k is the deterministic input to the battery (known); 
5 w k is the process noise (i.e. random inputs into the system); 

r k is the measurement noise (i.e. random noise that affects 
measurements); and 

A k , B k , F kj G k are the system matrices (known). 
[0023] The KF estimates a state by using a form of feedback 

10 control. In short, the KF estimates the state at some time and then 
obtains feedback in the form of noisy measurements. The KF 
equations fall into two groups: time update equations (predictor) and 
measurement equations (corrector). The time update equations project 
the current state and error covariance estimates to obtain estimates for 

15 the next time step. The measurement update equations incorporate 
new measurements into the estimates to obtain an improved estimate. 
[0024] The KF-based state estimator of the present invention 
introduces a mechanism that propagates the state estimate based on 
battery dynamics and propagates the estimate error covariance based 

20 on battery dynamics and knowledge of how noise affects the battery. 
The battery is a non-linear system. However, the KF-based state 
estimator simplifies the system by either linearizing non-linearities or 
removing them. Further, the KF-based state estimator accounts for 
time varying parameters to more quickly improve estimate errors. The 

25 controller processes the KF-based state estimator and controls other 
components of the system in response to the state estimates. 
[0025] The state estimates include SOC and internal (ohmic) 
resistance (R 0 hmic). The KF-based state estimator propagates 
(predicts) the estimated state value (x k ) and the covariance of the 

30 estimate error forward in time to the next measurement. The next 
measurement is obtained and the KF-based state estimator updates 
(corrects) X|< based on the measurement. The covariance of the 
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estimate error (P) is updated based on the measurement. As 
discussed in detail below, the order of these operations can vary. 
[0026] The process model of the KF-based state estimator 
describes how the internal states of the battery 12 change in response 
to current. The process model is provided as: 
x k+1 = Ax k + Bu' k + w k 



where: 



■[:;] 



(identity matrix); and 



10 B = 



Vce^Qsoc^i) 0 -1 
0 0 0 



1 (establishes relationship 
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between SOC and current using charge efficiency (Hce))- Expanding 
this vector equation, the following V 0 c and Rohmic process models are 
provided: 

[Voc]k + i = [V 0 c]k + Bu^ + w k 

1 5 [Rohmic]k+1 = [Rohmkjk + Bu' k + W k 

[0027] The measurement model of the KF-based state estimator 
describes how the states and the input current appear as a voltage 
across the terminals of the battery 12. The measurement model is 
provided as: 
20 y' k = F k x k + r k 

where: 

~ [l I scaled, k \ 

Expanding this vector equation, the following measurement model is 
provided: 

25 y' k = Voc + RohmicUcaled.k + T k 

[0028] At time k, the terminal voltage (V t ) and terminal current (l t ) of 
the battery 12 are measured using the meter 20. These 
measurements are provided as: 

yk = v t|k 
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u k = lt,k 

[0029] To include the effects of dynamics resulting from V d! and V hys t 
auxiliary relationships are implemented. Implementation of the 
auxiliary relationships enables the KF-based algorithm to be of a much 
5 lower order than would be typical. In this way, a more simple state 
estimator is provided. The dynamics are accounted for based on V t , k 
and l t ,k according to the following auxiliary relationships: 

Vhysi,k+\ ~ f^hyst^k » ^t,k ' ? SD,k ) 
I SD = f^QsOCk »^*) 

[0030] Compensated voltage (y') and current (u') values are 
10 determined based on V t ,k and l t ,k and the dynamics according to the 
following: 

yk' = [V tlk - V 0 - V d |, k - V hy st,k] 
7 a 

^scaled 
. J sd . 

V 0 is the open-circuit voltage or battery voltage at 0% SOC. 
1 5 [0031 ] Applying the compensated current value (u' k ) to the V 0 c and 
Rohmic process models provides: 

[Voc]k + i = [Voc]k +[lt,kncE(Qsoc,lt) - Isd]/3600 + w k 

[Rohmic]k+1 = [Rohmic]k + W k 

[0032] Applying the compensated voltage value (y' k ) to the 
20 measurement model provides: 

V tf k - V 0 — Vdi.k - Vhyst,k = Voc + Rohmic Ucaled.k + *k 

where: 

r k is the sum of all of the noises in the system. 
System noises include those associated with quantitization and A/D 
25 converter accuracy, those associated with the auxiliary equations and 
those associated with modeling errors between the states and 
measurements. 



7 



[0033] The estimated state values (Qsoc, Rohmic) are propagated 
according to the following equation: 

Xkik^ = Ax k -i|k-i + Bu' k (Process Propagation Model) 
The notation k|k-1 indicates the estimate of the value at time k given 
5 the information available at time k-1 . The Process Propagation Model 
(PPM) is based on the Process Model (PM) above. The process noise 
(w k ) drops out because only propagation of the mean is desired. The 
process noise (Wk) is randomly generated with a mean value of zero. 
[0034] An initial value for the battery states (Qsoc, Rohmic) are pre- 

10 programmed in memory. Certain assumptions are required to 

determine the initial value. It is appreciated that these assumptions 
can vary based on the characteristics of battery. Exemplary 
assumptions include the battery SOC = 50% fro SOC initialization. If a 
55 amp-hr battery is assumed, the nominal battery charge is equal to 

15 27.5 amp-hrs. If it is assumed that the battery 12 has an equal 

probability of having any charge within 15 amp-hrs of the 27.5 amp-hrs, 
then the variance is roughly 20(amp-hrs) 2 . The nominal resistance is 
assumed to be 3mohms. If it is assumed that the resistance has an 
equal probability of having any resistance within 1 .5mohms of 3mohms 

20 then the variance is 0.8(mohms 2 ). The assumptions are selected to 
minimize the error between the actual states and the assumed states 
to enable quicker error reduction. 

[0035] Assumptions are also required for the measurement and 
process noises. With regard to the measurement noise, there are four 

25 sources of noise in the measurement model: 1) the terminal voltage 
measurement, 2) the terminal current measurement, 3) the auxiliary 
equations, and 4) the measurement equation. Ignoring auxiliary 
equation and measurement equation noises, the terminal 
measurement variance (R) is estimated as 5(mV) 2 /cell. With regard to 

30 the process noise, values are chosen to achieve a desired behavior in 
terms of convergence to steady state and steady state blending of 
process model vs. measurement model corrections. Thus, process 
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model noise values (Q) were chosen as 5e" 5 (amp-hr) 2 and 6e* ! 
(ohms) 2 /cell. This provides the following matrices: 

R = [0.005] 

Tsio- 5 0 

0 610' 9 



5 It is appreciated that the initial values and assumptions described 
herein are merely exemplary in nature. The initial values and 
assumptions may vary based on the particular characteristics of the 
battery 12. 

[0036] A Kalman gain matrix (K) is calculated according to the 
1 0 following equation: 

Kk = Pkjk-iFk^FkPkik-iFk 1 + R)" 1 

where: 

F = [l I t k J (varies with time); and 

R is the measurement noise covariance (known). 
15 [0037] The covariance of the estimate error is propagated according 
to the following equation: 

Pkjk^ = APk-iik-iA* + Q (Error Propagation Model) 

where: 

A* is (known); and 
20 Q is the process noise covariance matrix (known). 

The Error Propagation Model (EPM) predicts the error covariance for 
the state estimate. 

[0038] The estimated state value (x k ) is corrected according to the 
following equation: 
25 Xk|k = Xk|k-i + Kk(y\ - FkX k |k-i) 

[0039] The covariance of the estimated error (P) is corrected 
according to the following equation: 

Pk|k = (Pk|k-1 1 + FkR" 1 Fk)" 1 

[0040] In the interest of reducing the amount of computing power 
30 and time required to process the KF-based state estimator, the KF- 
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based state estimator can be simplified. In particular, the A-matrix is 
eliminated to simplify the propagation of the state estimate and the 
propagation of the estimate error covariance. This provides the 
following simplified equations: 
5 Xkik-T = Ximiim + Bu' k (Simplified PPM) 

Pk|k-i = Pk-i|k-i + Q (Simplified EPM) 

The Kalman gain equation is simplified to eliminate the need for matrix 
10 inversion. This provides the following simplified equation: 

Kk = Pkjk-iFk^FkPkik-iFk 1 + R)" 1 

The error covariance update equation is also simplified, providing the 
following: 

15 Pk|k = (I - KkFk)Pk|k-1 

[0041] Referring now to Figure 4, the state estimator control will be 
described in detail. In step 100, control determines whether the system 
10 is on. If so, control continues with step 102. If not, control loops 

20 back. In step 102 control determines initial state values (xo|o) and initial 
error covariance values (Po|o). Selection of these values is described 
above. In step 104, control sets a counter, k equal to 1 . 
[0042] In step 106, control measures the terminal voltage and 
terminal current by receiving signals from the meter 20. In step 108, 

25 control determines the dynamic factors using the auxiliary 

relationships. In step 110, control determines the compensated input 
matrices. In step 112, control processes the KF-based state estimator 
to estimate the present SOC and resistance of the battery 12. Control 
processes the KF-based state estimator using the compensated input 

30 matrices. 

[0043] In step 114, control operates the system 10 based on the 
estimated SOC and resistance. Operation of the system 10 can 

10 



include battery charging, power control or system diagnostics. It is 
appreciated, however, that other operations can be included over those 
presently described. In step 116, control determines whether the 
system 10 is off. If not, control continues with step 118. If so, control 
5 ends. In step 118, control increments the counter by 1 then loops back 
to step 106. 

[0044] Referring now to Figure 5, the estimation control will be 
described in detail. The estimation control occurs within step 1 12 of 
Figure 4. In step 130, control predicts the current state estimates 

10 based on prior information (x k | k -i). In step 132, control predicts the 

current estimate error covariance based on prior information (P k | k -i). In 
step 134, control determines the current Kalman gain (K k ). In step 136, 
control updates the current state estimates (x k | k ) based on x k | k -i, Kk and 
y k \ In step 138, control updates the estimate error covariance (P k | k -i) to 

15 be used in the next iteration then control ends. 

[0045] The KF-based state estimator can operate using either 
backward rectangular integration, forward rectangular integration or 
trapezoidal integration. In order to design the KF-based state estimator 
to operate using backward or forward integration, the prediction and 

20 correction operations change order. For example, for backward 

rectangular integration, the prediction equations are processed first and 
the correction equations second. For forward rectangular integration, 
the correction equations are processed first and the prediction 
equations second. 

25 [0046] To enable trapezoidal integration the simplified PPM and the 
B-matrix are modified. In particular, the simplified PPM uses an 
average value of the compensated current values (u k ). This provides 
the following modified PPM equation: 

x k |k-i = Xk-i| k -i + B[(u' k + u' k )/2] 

30 The B-matrix is provided as: 
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B = 



WcE (QsOCk * A,* ) 



0 



0 
0 



-0.5 



VcE (QsOC,k-\ ' ?t,k-l ) 



3600 



[0047] The KF-based state estimator of the present invention 
provides the simplest state estimator for a dynamic system. The 
simplicity of the KF-based state estimator enables practical 
implementation in applications, such as automotive applications, 
because of the reduced amount of processing power and time 
required. The KF-based state estimator accounts for system dynamics 
and provides a linearized model of a non-linear system. 
[0048] Those skilled in the art can now appreciate from the 
foregoing description that the broad teachings of the present invention 
can be implemented in a variety of forms. Therefore, while this 
invention has been described in connection with particular examples 
thereof, the true scope of the invention should not be so limited since 
other modifications will become apparent to the skilled practitioner 
upon a study of the drawings, the specification and the following 
claims. 
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